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Abstract 

CO ■ In a previous paper we have shown that in quantum chromodynamics the gluon 

propagator vanishes in the infrared limit, while the ghost propagator is more 
singular than a simple pole. These results were obtained after angular averag- 
ing, but in the current paper we go beyond this approximation and perform an 
■ exact calculation of the angular integrals. The powers of the infrared behaviour 

of the propagators are changed substantially. We find the very intriguing result 
that the gluon propagator vanishes in the infrared exactly like p 2 , whilst the 
ghost propagator is exactly as singular as 1/p 4 - We also find that the value of 
the infrared fixed point of the QCD coupling is much decreased from the y-max 
pH | estimate: it is now equal to Att/3. 



^ Following a recent study by von Smekal et aLQ, we analyzed in Ref.Q the coupled Dyson- 

Schwinger equations for the gluon and ghost form factors F and G. The approximations were 
two-fold: firstly the vertices were taken bare, and secondly angular averaging was introduced 
(the so-called y-max approximation). Deferring to later work an improvement of the vertices, 
in this paper we seek to remove the deficiency of the y-max approximation. On the one hand 
the results might be regarded simply as quantitative adjustments to the y-max calculations; but 
on the other hand they are far from negligible. The numerical value of the infrared fixed point 
is reduced by a factor of almost three; and the finding that the gluon propagator has a simple 
zero, while the ghost propagator has a double pole, might perhaps be deemed a qualitatively 
new result. 

As an improvement on the y-max approximation used in Ref.|Q], we now solve the coupled 

integral equations for the gluon and ghost propagators with an exact treatment of the angular 
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integrals. Although the y-max approximation is good in the ultraviolet region, where the form 
factors run logarithmically, we will see that it is a crude approximation in the infrared region, 
where the form factors exhibit power behaviour. The y-max approximation apparently leads to 
a difficulty in the ghost equation, as the asymmetry in the treatment of the gluon and ghost 
momenta in the loop gives an ambiguous result. In our previous study we used the ghost 
momentum for the radial integration and the gluon momentum for the angular integration. 
However, if we exchange these momenta we no longer find a consistent solution to the integral 
equation. We will show later why the first choice is a better approximation than the latter, 
thus motivating the choice made in Ref.[||. Furthermore we will show that both choices are 
equivalent if we treat the angular integrals exactly, although this equivalence is by no means 
trivial. 

With the radial integration over the ghost momentum, the ghost form factor satisfies 

1 - z 3 - — [ A dy yG(y) f d9 sin 4 9 ^ , (1) 



G(x) TT Jo Jo 

whereas with the radial integration over the gluon momentum, we have instead 

1 -Js-- f^dyFiy) f^sin 4 ^, (2) 



G(x) TT Jo Jo Z 

where A = g 2 /IQtt 2 , g being the strong coupling constant and where z = x + y — 2^/xycosO. 
Here we have used the fact that Z\, the ghost-gluon vertex renormalization constant, is equal 
to unity in the Landau gauge. 

We will show that 

F(x) = Ax 2k G(x) = Bx- K (3) 

are solutions of Eq. ([l|) and Eq. (|2|), yielding identical conditions on XAB 2 , but that this equiv- 
alence does not hold in the y-max approximation. 

Substituting Eq. (|3]) in Eq. (|T|) and evaluating the angular integral in terms of the hypergeometric 
function, 2-F1, we obtain 



= Zs-^L^fdyy-^^ (-2« + 2,-2k;3;| 
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+ / dyy K - 1 2 F 1 I -2k + 2, -2k; 3; - 

x V y 



Set t = y/x in the infrared integration and t = x/y in the ultraviolet integration and take 
A 2 — > oo: 



x 
~B 



9XAB 



j- I dt ( t~ K+1 + t-^ 1 ) 2 Fi (-2k + 2, -2k; 3; t) . 



Note that we have dropped Z3. This is in accordance with standard regularization procedure: 
the integral in Eq. (|lj) is convergent if Re k < 0, whereas a subtraction is necessary if < k < 1. 
By identifying Z3 with this subtraction constant, we ensure that G(x) is defined by analytic 
continuation in k beyond Re k = 0. This continuation is made explicit in terms of the generalized 
hypergeometric function 3.F2 (see Eq. (|l~7|)). After matching the coefficients of x K we obtain 



1 



XAB 2 



1 



2 — K 
1 



]F 2 (-2k + 2, -2k, 2 - k; 3, 3 - k; 1) 



(5) 



F 2 (-2k + 2, -2k, -re; 3, 1 - k; 1) 



The y-max approximation that we made in [Q] amounts to replacing /c by zero in the hypergeo- 
metric function in Eq. (H). The result was 
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XAB 2 



1 1 



2 — K K 



(6) 



Next we make a similar analysis of Eq. (|2|): 

9XAB f A2 
4 Jo 



x 



dyy 2K y^-\F 1 (K+l,K-l-?,; V 



(7) 



B 4 jo \ y>, 

After defining t = y/x in the infrared integration and t = x/y in the ultraviolet integration, 
taking A 2 — > 00 and eliminating Z% as before, we obtain 
x K 9AAB 



B 



, -x 

4 jo 



K (t 2K + t""" 1 ) 2 Fi (k + 1, /c - 1; 3; t) . 



(8) 



Using Eq. (pL7|) in Eq. (0) and equating the coefficients of x K we find 



1 



AA£ 2 



1 



3 F 2 (k + 1, K - 1, 2k + 1; 3, 2k + 2; 1) 



2k + 1 

3F2 (k + 1,k - 1, -k; 3,1 — k; 1) 

K 



(9) 



The y-max approximation of this expression is once more obtained by replacing k by zero in the 
hypergeometric function in Eq. (H), yielding 



1 
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XAB 2 



1 1 

+ 



2k + 1 3(2k + 2) k 30-1) 



(10) 



This approximation is quite different from that of Eq. (^) , illustrating the precariousness of the 
y-max approximation. 

However, the exact forms Eq. (|J) and Eq. (^) are identical, as we verified by using Mathematica 
to check the numerical equivalence over a domain of k values. 

We now consider the gluon equation, and will show that it too is solved by a pure power behaviour 
if we keep only the ghost loop. Even in the presence of the gluon loop, when the power behaviour 
is no longer an exact solution of the equation, it will still represent the correct leading order 
infrared asymptotic behaviour. Omitting then the gluon loop, we write the equation for the 
gluon form factor as follows: 

= Z 3 + - -JL G (y) / d9 sin 2 6 M(x, y, z)G(z) , 
Fyx) TT Jo X Jo 

where 

/ ' x + y y 2 \ i i 2 y z 

M(x,y,z)= — \- + - + JL --. 

\ 2 x z 2 x x 



id 



We now substitute the solution (||) into Eq. (11), obtaining 

2\ 



X 



- 2 " 2XB 2 f A2 dy 



A 



IT JO X 



d9 sin z 6 



x + y y* 



-K-l 



+ l 2 + l' ■■" 



We evaluate the angular integral with the help of Eq. (16): 



x 



-2k 



A 



Z 3 + XB 2 / A %-* 
Jo x 



x + y _ r_ 

2 x 



y> K - 1 2F 1 (« + 1,k;2; 



y< 
y> 



+ i] + f)^(-- 1 ^ 



] ^— 2 F 1 f«-l,«_2;2;^ 
x \ y> 



We split the integration region and define t = y/x in the infrared integration and t = x/y in the 
ultraviolet integration, as before, and finally implicitly implement the analytic regularization, 
taking A 2 — > oo. After some rearrangement we obtain 



^— = XB 2 x- 2k I dt 
A Jo 



(12) 



| [i (f- K + t~ K+1 + t 2K " 1 + t 2K ~ 2 



^— k+2 _|_ 3 



+ 



1 (+-k + t 2K-2^j + 2 A-ie+1 + t 2 K -3 



2 F 1 (K + l,K;2;t) 
2 F 1 {n,n-\-2-t) 



t~ K + t 2K ~ 3 



2 F 1 (K-l,K-2;2;t)\ . 



4 



Using Eq. p7| ) in Eq. ( pL2[ ) and equating coefficients of x K , we find 
1 



(13) 



XAB 2 

2(1 ] _ k) 3^2 (« + 1, «, 1 - is; 2, 2 - «; 1) + T^^y 3^2 (k + 1, «, 2 - k; 2, 3 - k; 1) 

+-/- 3^2 (« + 1, k, 2k; 2, 2k + 1; 1) + — 1 3^2 (k + 1, k, 2k - 1; 2, 2k; 1) 
4k 2(2k — 1) 

--— 3 F 2 (#c + 1, k, 3 - k; 2, 4 - «; 1) - 1 3^2 (« + 1, «, 2k - 2; 2, 2k - 1; 1) 
6 — K 2{K — 1) 

+ wrr \ 3 F 2 (k, k - 1, 1 - k; 2, 2 - k; 1) + 1 3^2 (k, k - 1, 2k - 1; 2, 2k; 1) 
2(1 — kJ 2(2k — 1) 

+ 77^— 3^2 (k, k - 1, 2 - k; 2, 3 - k; 1) + — - 3 F 2 (k, k - 1, 2k - 2; 2, 2k - 1; 1) 

2 — K K — 1 

-— — 3^2 (k - 1, k - 2, 1 - k; 2, 2 - k; 1) - - 1 3^2 (k - 1, k - 2, 2k - 2; 2, 2k - 1; 1) . 

1 — K 2(k — 1J 

Some of these generalized hypergeometric functions are in fact infinite for certain values of k; 
they should be understood by replacing the final argument '1' in each generalized hypergeometric 
function by r, and then by taking the limit r — ► 1. Since this amounts to replacing the upper 



integration limit in Eq. (12) by r, the procedure is clearly legitimate. 

Finally, to demonstrate that Eq. (^) really is a solution of the coupled integral equations, for a 
particular value of k, we equate the right-hand side of Eq. (||) or Eq. (||) with that of Eq. ( Jl3| ) 
and solve for k. We used the routine FindRoot of Mathematica and discovered that, in the limit 

K — ► 1 , 

both expressions become equal. This limit has been checked analytically in App. |B|. We find 

XAB 2 = - 
3 

which means that the infrared fixed point for the gauge invariant running coupling, a(p 2 ) = 
An\Z 2 F(p 2 )G 2 {p 2 ), has the value 

a(0) = y « 4.19 

which is quite different from the the value we found using the y-max approximation. There we 
had a(0) ~ 11.47, nearly three times the exact result, corresponding to k ~ 0.77. 

Our findings for the propagators in the infrared region can be summarized by the following 
formulas: 
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duon propagator in Landau gauge 



<(P) ~ -5° 



P 



ghost propagator in Landau gauge 



G ab (p) ~ -5 ab — 



A Angular and Radial Integrations 



Throughout we used the angular integration formula 



dd sin 2r 6 z n = B ( r + -, - ] 2 F± ( -n, -n - r; r + 1; — 
o V 2 2/ V y > 



The radial integrations of hypergeometric functions are evaluated as follows: 

dtt 



/ dtt v 2 F 1 (a,b;c;t) = T 
J o „ =0 

oo 

= E 



(a) n (b) n f 1 ^ + v + n 



(a)n(b) n 1 

. (| (c) n n! v + n+l 



1 



^ + 1 



3 F 2 (o,6,i/ + 1;c,i/ + 2;1) 



B Limit of XAB 2 for « -> 1 



-~{ 3 F 2 (0, -2, 1; 3, 2; 1) - Km 3^2 (2 - 2k, -2, -1; 3, 1 - «; 1) J 



We check analytically that the limit of XAB 2 for k — > 1 is identical for Eqs. ( |l3| , 
from Eq. (||) we find 

1 

XAB 2 
and from Eq. (||) 

= -^{^ 3^2 (2, 0,3; 3, 4; 1) - lim 3 F 2 (2, k - 1, -1; 3, 1 - «; 1) 



1-1 + 
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We make an analoguous check for the gluon equation. To circumvent convergence problems, we 
replace the argument of the generalized hyper geometric functions by r, and take the limit t — ► 1 
at the end. From Eq. (|T3|), 

1 

XAB 2 ~ 

lim lim ( 1 3^2 (2, 1, 1 - «; 2, 1; r) + 3 F 2 (2, 1, 1; 2, 2; r) - - 3 F 2 (2, 1, 2; 2, 3; r) 

r^l K^l 1.2(1 — K) 4 

— 7T7 — ~ — 7\ 3^2 (2, 1, 2k - 2; 2, 1; r) + — 1 3^2 (1, k - 1, 1 - «; 2, 1; r) 
z(k — 1J 2(1 — Kj 

+~ 3^2 (1, 0, 1; 2, 2; r) + — !— 3 F 2 (1, k - 1, 2k - 2; 2, 1; r) 
2 re — 1 

-— !— 3^2 (« - 1, -1, 1 - «; 2, 1; r) - 1 3^2 (« - 1, -1, 2k - 2; 2, 1; r) } 

1 — K 2(K — Ij J 

, , f 1 (l-K)n + (2K-2) n _ , ,1 , . I 5 

= £3 { 2(T^o S 1 -- n + ^ (i. 1; 2 =-) - 4 2Fl & 2 ^} + 2 

y ' n=l J 

I ~J V 2n + n + 1 2(2 + n) J T J + 4 
t— 1 l^A 2n n 2n/ J 

= 3. 
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